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1 Introduction 


Charged Higgs bosons go along with many extensions of the Standard Model, such as supersymmet¬ 
ric versions of the Standard Model or general Two-Higgs-Doublet models. The neutral Higgs-like 
particle with a mass ~ 125 GeV, discovered by the ATLAS and CMS experiments [1, 2], behaves 
within the presently still sizeable experimental uncertainties like the Higgs boson of the Standard 
Model (see [3, 4] for latest results), but on the other hand leaves ample room for interpretations 
within extended models with a richer spectrum. A scenario of particular interest thereby is the 
Minimal Supersymmetric Standard Model (MSSM) with two scalar doublets accommodating five 
physical Higgs bosons, at lowest order given by the light and heavy CP-even h and H, the CP- 
odd A, and the charged Higgs bosons. The discovery of a charged Higgs boson would constitute 
an unambiguous sign of physics beyond the Standard Model, providing hence a strong motivation 
for searches for the charged Higgs boson. 

Experimental searches for the charged Higgs bosons of the MSSM (or a more general Two-Higgs- 
Doublet Model) have been performed at LEP [5] yielding a robust bound of > 80 GeV [6] . The 
Tevatron bounds [7] are meanwhile superseeded by the constraints from the searches for charged 
Higgs bosons at the LHG [8]. 

The Higgs sector of the MSSM can be parametrized at lowest order in terms of the gauge 
couplings gi and 52 , the mass itia of the CP-odd Higgs boson, and the ratio of the two vacuum 
expectation values, tan/3 = V 2 /V 1 ; all other masses and mixing angles are predicted in terms of 
these quantities. Higher-order contributions, however, give in general substantial corrections to the 
tree-level relations. 

The status of higher-order corrections to the masses and mixing angles in the neutral Higgs sector 
is quite advanced. A remarkable amount of work has been done for higher-order calculations of the 
mass spectrum, for real SUSY parameters [9-23] as well as for complex parameters [24-30]. They 
are based on full one-loop calculations improved by higher-order contributions to the leading terms 
from the Yukawa sector involving the large top and bottom Yukawa couplings. Quite recently, 
the C(at) terms for the complex version of the MSSM were computed [29, 30]; they are being 
implemented into the program FeynHiggs [31-33]. 

Also the mass of the charged Higgs boson is affected by higher-order corrections when expressed 
in terms of itia- The status is, however, somewhat less advanced as compared to the neutral Higgs 
bosons. Approximate one-loop corrections were already derived in [34-36]. The first complete 
one-loop calculation in the Feynman-diagrammatic approach was done in [37], and more recently 
the corrections were re-evaluated in [28, 38, 39]. At the two-loop level, important ingredients for 
the leading corrections are the 0{atas) and 0{a1) contributions to the charged self-energy. 
The 0{atOis) part was obtained in [27] for the complex MSSM, where it is required for predicting 
the neutral Higgs-boson spectrum in the presence of CP-violating mixing of all three neutral CP 
eigenstates with the charged Higgs-boson mass used as an independent (on-shell) input parameter 
instead of niA- In the CP-conserving case, on the other hand, with ma conventionally chosen 
as independent input quantity, the corresponding self-energy contribution has been exploited for 
obtaining corrections of 0{atas) to the mass of the charged Higgs boson [39]. In an analogous 
way, the recently calculated 0{a^) part of the self-energy in the complex MSSM [29, 30], can 
now be utilized for the real, CP-conserving, case to derive the 0{al) corrections to the charged 
Higgs-boson mass as well. 

In the present paper we combine the new two-loop terms of 0{al) with the complete one-loop 
and O(atas) two-loop contributions to obtain an improved prediction for the mass of the charged 
Higgs boson. The results have been implemented into the code FeynHiggs. An overview of the 
calculation is given in Section 2, followed by a numerical evaluation and discussion of the two-loop 
corrections in Section 3 and Gonclusions in Section 4. 
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2 Higgs-boson mass correlations 

2.1 Tree-level relations 

We consider the Higgs potential of the MSSM with real parameters, at the tree-level given by 

fdiiggs = + 1712 'H\T-L2 + C-) 

+ ^ (5? + 9 l ) (^2^2 - hIhi)' + iff! (nln2) {ulUi 

with the mass parameters m|, 771^21 the gauge-coupling constants gi,g 2 - The two scalar 
Higgs doublets in the real MSSM can be decomposed according to 
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with real vacuum expectation values vi and V 2 - The ratio U 2 /U 1 is denoted as tan/? = tp. The 
mass-eigenstate basis is obtained by the transformations 
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[with Sx = sin x and Cx = cos x ], where h, H, A and denote the physical neutral and charged 
Higgs bosons, and G°,G^ the unphysical neutral and charged (would-be) Goldstone bosons. 

The Higgs potential in the real MSSM can be written as the following expansion in terms of the 
components h, H, A, ,G^ [with (i?“)^ = (G“)'^ = G+], 
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omitting higher powers in the field components. Explicit expressions for the entries in the mass 
matrices are given in Ref. [28] for the general complex MSSM [the special case here is obtained for 
setting Ta = 0 in those expressions]. Of particular interest for the correlation between the neutral 
GP-odd and the charged Higgs-boson masses are the entries for m\ and , reading 

m\ = ml + ml + m ^2 S2/3„ - 3 ( 5 ? + gDivf - vj) C 2 p„ , 

= ml s\ + ml + ml2 S2p, - \{gl + gDivf - vl) C2p, + \gl(viCp^ + V2Sp^f . 

At lowest order, after applying the minimization conditions for the Higgs potential, the tadpole 
coefficients Th^Tn vanish and the mass matrices become diagonal for Pc = Pn = P, yielding 

TO^± =m\ + , 
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when a is chosen according to 
tan(2a) = 


m\ 


■ m^ 


m\-mz 

The Goldstone bosons G° and G^ remain massless. 

In the following we focus on the modification of the relation (2.6) by higher-order contributions, 
which allows to derive the charged Higgs-boson mass in terms of the A-boson mass and the 
model parameters entering through quantum loops. 
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2.2 The charged Higgs-boson mass beyond lowest order 

Beyond the lowest order, the entries of the mass matrix of the charged Higgs bosons are shifted 
by adding their corresponding renormalized self-energies. The higher-order corrected mass M[j± 
of the physical charged Higgs bosons, the pole mass, is obtained from the zero of the renormalized 
two-point vertex function. 


Mjj± — 3fie(so), ^H+H- (p^) 


= i 


V =.50 


/ - TO^± -f (p^) 


= 0 . 


V =-50 


(2.9) 


Therein, "Eh+h- (p^) denotes the renormalized self-energy for the charged Higgs bosons which 
we treat as a perturbative expansion. 
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At each loop order fc, the renormalized self-energy is composed of the unrenormalized 

self-energy and a corresponding counterterm according to 

(P‘) = (p>) - (pP ). (2.11) 

At the one-loop level the counterterm is given by 

(p^) = (m^± - p^) , (2.12) 


and at the two-loop level by 


= (m^±-p^) d^‘^'^ZH±H± + I (5 ^^^Zh±h±) 
+ 5^^'^Zh±h± + i 5^^^Zh±g± 


-p2i (2.13) 


involving field-renormalization constants and genuine mass counterterms of one- and two-loop order; 
they are specified in Ref. [30], from where conventions and notations have been taken over and 
simplified to the case of the real MSSM. 

Whereas the one-loop self-energy S^+^_(p^) of the charged Higgs boson is completely known, 
at the two-loop level only results in the approximation for = 0 have become available, namely 
the 0{atas) corrections calculated earlier [27, 39], and the two-loop Yukawa contributions 0(aj) 
which are presented in this paper. The evaluation of these terms is performed in the gaugeless limit 
and the bottom-quark mass set to zero (as done in Ref. [39]), thus yielding the top-Yukawa-coupling 
enhanced parts. Detailed analytical results of the two-loop self-energy and renormalization were 
published in Ref. [30]. The diagrammatic calculation of the self-energies and counterterms was 
performed with FeynArts [40], FormCalc [41], and TwoCalc [42]. The full list of Feynman diagrams 
of 0(a1) for the self-energy of the charged Higgs boson is illustrated in Fig. 1. 

Within our approximations for the two-loop part of the charged Higgs-boson self-energy, 

(0) = (0) - (0), (2.14) 


the two-loop counterterm (2.13) simplifes to 
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Figure 1. Full list of two-loop self-energy diagrams for the charged Higgs bosons. Each cross denotes a 
one-loop counterterm insertion. = h, H, A, G; = H~, G~. 


The genuine mass counterterms are determined by Eq. (2.5) and setting /3„ = = /3 (see 

also Ref. [30]). In the gaugeless limit they are given by (for k = 1,2) 

. (2.16) 

The other genuine mass counterterms are determined by the relation 

= ~ 2s ~ 

involving the tadpole counterterm and the counterterm for the renormalization of 

tan /3. 

In the real MSSM, the mass of the CP-odd Higgs boson niA is conventionally chosen as a free 
input parameter; it can thus be renormalized on-shell at each order. Accordingly, the corresponding 
renormalization conditions in our present approximation read in terms of the renormalized A-boson 
self-energy as follows, 

E^*^^(0) = E^^^(0)-5('=)m5(0) =0. (2.18) 
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Figure 2. Full list of two-loop self-energy diagrams for the A-boson. Each cross denotes a one-loop 
counterterm insertion. = h, H, A, G; "hp = A, G; = H~, G~. 


( 2 ) 

The unrenormalized self-energy corresponds to the Feynman diagrams depicted in Fig. 2. The 
counterterms in (2.18) at the one-loop and two-loop level read as follows, 


= m\ 5 ^^'^Zaa + ■ 




AA 




+ S^’^^tua 


^AG ■ 


S^^^mA(O) = m\ 

+ 5^^^Zaa + 5^^^Zag 

The one-loop non-diagonal mass counterterm therein is given by 

S^^'^Th - m\ c| . 




2 Sw Mw 


(2.19a) 
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( 2 . 20 ) 


From the conditions (2.18) for fc = 1,2 the renormalization constants are determined and 

thus the mass counterterms for the charged Higgs bosons in Eq. (2.16), required for 

the two-loop counterm (2.15) in the charged Higgs-boson self energy. All field-renormalization con¬ 
stants S^^'^Z^aa,ag,h^h±.h^G^} s-re linear combinations of the basic field-renormalization constants 
6^^'>Z'fi. for the two scalar doublets (2.2), as given in Ref. [30]. 

In addition to the mass counterterms the independent renormalization constants re¬ 

quired for renormalization of the charged Higgs-boson self-energy are: the field renormalization 
constants S^^'^Z-f^., the renormalization constant for tan/3, and the tadpole renormalization 

constants S^^'^Th, S^^'^Th (the two-loop field renormalization constants cancel in the renormalized 
self-energies in the = 0 approximation). Moreover, for the one-loop subrenormalization, we need 
the counterterms for the top quark and squark masses ^(1) 

mj.^, and for the trilinear 

coupling as well as the counterterm for the bilinear coefficient of the superpotential, 

They are fixed in the same way as described in Ref. [30] and we do not repeat them here. 
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3 Numerical analysis 

In this section we compute numerically the charged Higgs-boson mass Mfj± in the real MSSM in 
terms of niA chosen as an input parameter. For this purpose, we combine in the renormalized 
charged Higgs-boson self-energy our new contribution described in the previous section with 

the already known complete one-loop term and the 0{atas) contribution, 

(/) = (/) + (0) + (0), (3.1) 

as the currently best approximation for (2.10). The resulting charged Higgs-boson mass Mh± is 
obtained via Eq. (2.9) with the help of FeynHiggs. 

In the following numerical analysis we use the input parameters as listed in Tab. 1 (giving also 
those parameters not needed for the two-loop self-energies, but required for specifiying the input 
for the other terms in (3.1) and for FeynHiggs). The other parameters of the MSSM not contained 
in Table 1 are kept variable and are given in the figures. The quantities /i, tp and the Higgs field- 
renormalization constants are defined in the DR scheme at the scale mt (see also Ref. [30] for more 
details). 


Table 1. Default input values of the MSSM and SM parameters. 


MSSM input 

SM input 

M2 

= 200 GeV, 

mt = 173.2 GeV, 

Ml 

= (54)/(3Cw) ^# 2 , 

mb = 4.2 GeV, 

"iq = 

= 2000 GeV, 

mr = 1.777 GeV, 

W 91 = rUu^ = 1713 

= 2000 GeV, 

Mw = 80.385 GeV, 

II 

II 

3 

= 0 

Mz = 91.1876 GeV, 


= 2000 GeV, 

Gf = 1.16639- 10-^ 

niq^ = rric^ = 

= 2000 GeV, 

= 0.118, 

a. 

II 

II 

0 

= 0 . 



The influence of the O (a^) corrections on the charged Higgs-boson mass decreases with increasing 
values of tp, where the top Yukawa coupling is diminished. Therefore we constrain our analysis on 
values of tp < 10. In the case of larger tp also the corrections of 0{ai,at) may become relevant (see 
also Ref. [39] for more discussions on the validity range). 

The shifts in the charged Higgs-boson mass resulting from the 0(aj) contributions are in general 
small. In Fig. 3 the dependence of M[j± on the Higgs-sector input parameter mA and on the third- 
generation soft-breaking squark mass parameter = rriq^ = rrij:^ = is depicted, showing a 
decreasing size of the two-loop mass shift (red) for increasing values of both variables. The upper 
section of the figure shows the charged Higgs-boson mass as obtained at the one-loop level (dashed), 
and with the inclusion of the 0(atas) contributions (green) and also the terms (blue). The 

lower section of Fig. 3 shows the mass shift originating solely from the 0{af) two-loop part. 
Thereby, the ©(a^) corrections appear as negative, thus diminishing the two-loop contribution 
of 0{atas). In total, the two-loop terms still yield a positive shift upon the one-loop result for Mfj±. 

Fig. 4 contains the charged Higgs-boson mass Mh±, together with the two-loop shift of 
for a typical low-mn scenario (left) [43] and for a scenario with heavier (right), versus the 
Higgsino mass /x. For large values of p,, the charged Higgs-boson mass Mh± decreases, but the 
mass shift AMh± resulting from the 0[a^) contributions becomes more sizeable, reaching I GeV 
and more for the low M^± case. In the scenario shown in the right panel of Fig. 4 the two-loop 
contributions are smaller in comparison to the one in the left panel, which is a consequence of the 
smaller Yukawa couplings for larger values of rriA and tan/3. 
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— without 0{aj), 0(att 
m~f = 1000 GeV 
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rriA [GeV] 


_ without 0(aj) 

;) _ mass shift by 0(aj) 
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Figure 3. Upper parts: prediction for the charged Higgs-boson mass M[j± including all known contribu¬ 
tions (blue), without the ©(q?) contributions (green) and without any two-loop corrections (black dashed) 
[mjj± is the tree-level mass according to Eq. (2.6)]. Lower parts: the mass shift AMjj± by the 0[at) con¬ 
tributions (red). Left: ruf = rriq^ = = 1000 GeV. Right: itia = 500 GeV. The other input 

parameters are tp = 8, ^ = 2000 GeV, = mf^ = 1000 GeV, Xt = 2 m^, Ab = At — 0, rrig = 1500 GeV, 
for both cases. 


tl3 - 7 , 


— total — without 0(aj) 

— without 0{aj), C)(Qr,Q;s) — mass shift by 0{a^) 
ma = 120 GeV 




ju [GeV] 


500 1000 1500 2000 2500 3000 3500 4000 
ju [GeV] 


Figure 4. Upper parts: prediction for the charged Higgs-boson mass M[j± including all known contribu¬ 
tions (blue), without the ©(a?) contributions (green) and without any two-loop corrections (black dashed). 
Lower parts: the mass shift AMu± by the ©(q?) contributions (red). Left: tp = 7, tua = 120 GeV, 
At = 2.5 rriq^, At = Ar = 0. Right: tp = 8, mriA = 500 GeV, Xt = 2mq^, Af, = At = 0. The other input 
parameters are rriq^ = rni^ = = 1000 GeV, = 1000 GeV, rrig — 1500 GeV, for both cases. 
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m-f 

Figure 5. Mass Mu± of the charged Higgs boson with all available higher-order terms and without the 
0[at) contributions, for ranges of mj and Xt allowed by HiggsBounds and mh = 125 GeV. Other parameters 
are niA = 200 GeV, tp = 8, fi = 3000 GeV, Ah = Ar = 0, = nif^ = 1000 GeV, nig — 1500 GeV. 

In all cases, the 0(^aj) contributions appear with negative sign and reduce slightly the positive 
mass shift arising from 0{atas). In general, the combined two-loop corrections result in a positive 
shift, which can amount to several GeV, on top of the one-loop prediction for M^±. 

In the figures mentioned above, the constraint rrih = 125 ± I GeV on the light Higgs-boson mass 
is imposed, except for the Iow-toj^ scenario in Fig. 4 (left) where it is the heavier 7J-boson that 
appears with a mass around 125 GeV (a scenario which may soon be excluded by more stringent 
limits on the charged Higgs-boson mass). One has to keep in mind, however, that not all of the 
parameter values in the figures, which are shown for illustrating the parameter dependence, will 
actually be allowed when more comprehensive phenomenological studies on the properties of the 
Higgs particle at 125 GeV will be performed. We have added such a more comprehensive analysis 
by probing the regions compatible with the experimental constraints by means of the program 
HiggsBounds [44-46]. The result is shown in Fig. 5, where the 0(^a^) effects for are displayed 
for possible combinations of the stop-sector parameters. Also here we find negative mass shifts in 
the typical range from —0.5 GeV to —0.8 GeV. 

4 Conclusions 

We have calculated the two-loop 0(^a^) contributions to the mass of the charged Higgs 

boson when derived from the A-boson mass ttia as an on-shell input parameter within the real, 
CP-conserving, MSSM and combined them with the complete one-loop and the two-loop 0{al') 
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contributions. We have presented numerical studies for scenarios of current phenomenological 
interest and discussed the effects of the various two-loop terms. 

The 0(^af^ two-loop corrections appear with opposite sign and smaller size with respect to the 
0{atas) contributions; in combination, the two-loop terms yield a positive shift to the mass of the 
charged Higgs boson as calculated at one-loop order. This shift in can be at the level of 

several GeV and thus of a size that may be relevant for the LHC (and a future electron-positron 
collider). 

The set of two-loop corrections considered here are expected to be particularly relevant in pa¬ 
rameter ranges of the real MSSM where the top-Yukawa terms provide a good approximation to 
the complete one-loop result, especially for relatively low values of tan/3 and tua- In this range, 
besides precise mass predictions, the experimental constraints on the mass and the phenomenolog¬ 
ical features of the lightest Higgs are important and play a substantial role when comprehensive 
analyses withing the MSSM Higgs sector are performed. 

Our results for the charged Higgs-boson mass have become part of the Fortran code FeynHiggs. 
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